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Abstract. Asymptotic Giant Branch (AGB) stars are important players in 
the chemical evolution modelling of galaxies, because they are major producers 
of several chemical elements and excellent tracers of the structure and of the star 
formation activity of their parent galaxies. A few examples on the importance of 
AGB stars are presented in this review, together with a number of open problems 
affecting chemical evolution model predictions related to the element enrichment 
C*") . by AGB stars: the evolution of 4 He and Na in globular clusters, the evolution 

of 3 He and the carbon isotopes in the Galactic disk, and the evolution of N and 
O in different types of galaxies. The need of homogeneous and complete sets of 
yields is emphasized. 
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<^ ■ 1. Introduction 

23: 

Qh' There are several reasons why galaxies and galaxy modellers care about Asymp- 

totic Giant Branch stars. AGB stars trace the halo stellar distribution (see 
Demers in this volume), thus indicating what is the size and the dynamics of 
the system; they trace intermediate-age stars (see Grebel, this volume), thus 
showing if and how much star formation (SF) activity has occurred in the re- 
gion 0.5 - 1 Gyr ago; they are the major site of production of several chemical 
elements (see Busso and Lattanzio in this volume) and fundamental contributors 
to the stellar nucleosynthesis yields, which, in turn, are the main ingredients of 
chemical evolution models. 

This review is focussed on the third point, which is the most important one 
for chemo-dynamical models of galaxy evolution. Let me however mention how 
important it is to trace intermediate age and old populations. IZwl8 is the most 
metal poor star-forming galaxy ever discovered. Being very blue and full of gas, 
it was often considered a local counterpart of primeval galaxies and a really young 
system, with SF activity begun only a few Myr ago. Ho wever, when HST optical 



phot ometry became available, AGB stars were found (jAloisi. Tosi. fc Greggio 
1999) and then confirmed with near-infrared HST photometry ( Ostlin 2000). 
The presence of AGB stars undisputably demonstrated that IZwl8 is not as 
young as originally thought and must have started forming stars at least 0.5-1 
Gyr ago. 

Here, I will describe the effects of AGB stars as interstellar medium (ISM) 
polluters with reference to the evolution of a few interesting elements: 3 He, 4 He, 
12 C, 13 C, N and Na. The galaxy evolution models described here are all standard, 
in the sense that they do not explicitely treat dynamical aspects. It is however 
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becoming increasingly clear that to understand galaxy evolution dynamical and 
hydrodynamical effects cannot be left aside. 



2. AGB stars as ISM polluters: 4 He, Na and the evolution of glob- 
ular clusters 

In recent years there has been an increasing interest on the possibility of a 
second generation of stars in some globular clusters. Some of the observed prop- 
erties of globular clusters are indeed considered evidence of a second SF event. 
For instance, the anticorrelation betwee n the sodium and oxygen abundances 
measured in stars of several clusters (e.g. ICarretta. Bragaglia. Cacciari 1120041 . 



and references therein) has been interpreted as the consequence of the cluster 
self-enrichment, if the gas replenishment from the retention of stellar ejecta is 
sufficient to allow for further SF. In this scenario, the stars born during the sec- 
ond SF episode form from mostly (if not entirely) recycled gas and their initial 
chemical compositions r eflect the yield s of th e star s mostly contrib u ting t o the 
cluster self-enrichment. iGratton et al.1 ( 200ll ) and iD'Antona et al.l ( 20021 ) sug- 



gested that high Na enrichment and significant O depletion are most naturally 
explained if the major culprits of the cluster self-pollution are relatively mas- 
sive AGB stars. The stars with more O and less Na would be those formed in 
the first generation, while the st ars with less Q and m ore Na would be those 
formed in the second generation. iDAntona et al. ( 20021 ) further suggested that 



the second generation would be significantly enriched also in helium (again as 
a consequence of the predominance of AGB star ejecta in the gas available for 
SF) and that this could explain the Horizontal Branch morphology of clusters 
with extreme blue tails. Other authors, however, have argued that the winds 
of massive stars have more chances than AGB stars to adequately pollute the 
cluster's medium without the sid e effects of requiring unusual in itial mass func- 
tions and stellar remnants (e.g. iPrantzos Charbonnell 120061 ). Appropriate 
chemo-dynamical models are required to test advantages and disadvantages of 
the various hypotheses. 

The discovery of a second, bluer Main Sequence (MS) in uj Centauri (jBedin et al 



2004), together with that of multiple subgiant and red giant branch (RGB) se- 



quences, also calls for the existence of multiple stellar generations, with the 
additional striking surprise, provided by high-resolution spectrosc opy, that the 



bluer MS is 0.3 dex less metal poor than the standard red MS (jPiotto et al 



20051 ). Piotto et al. argue that the only way to allow for the measured colour 



shift between the two MSs with the measured metallicity difference is to let the 
bluer MS be much more helium rich than the other, with a difference in the 
helium mass fraction AY=0.14. If we attribute to the red MS a primordial He 
mass fraction of Y=0.24, this implies that the blue MS should have Y=0.38: an 
abundance higher than in any other observed star cluster or galaxy ! 

Would AGB stars be able to provide such a huge helium enhancement ? It 
is very unlikely. However, uj Centauri is definitely not a normal globular cluster; 
may be not a cluster at all, but the remnant of a nucleated dwarf galaxy captured 
and stripped by the Milky W ay, with the current c luster actually being the 
original nucleus of the satellite. Bekki <fe Norris I ( 200(1) have recently shown that 



the observed properties of u Centauri cannot be explained without considering 
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the strong dynamical interactions with the Galaxy. The actual question then is 
whether the medium, out of which subsequent stellar generations have formed, 
was enriched by u Centauri own stars or by the stars of the host dwarf originally 
surrounding it. Whether the polluters were AGBs, massive stars or supernovae 
is in this case a second level issue. 
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Figure 1. Predictions from chemical evolution models (solid lines) for 
w Centauri compared with the corresponding observational data (dots 
and dashed l i nes). Left panel: the stellar metallicity distribution from 
Solli ma et alJ (120051) : central top panel: age-metallicity r elation from 
Hilke r et alj (|2004f ): central bottom panel: helium vs iron from iPiotto et al.l 
(2005|, box) and in RR Lyraes (dots; from lSollima et ai1l2006f ): right panel: 
abundance ratios vs iron from several literature sources. 



We are currently modelling the chemical evolution of uj Centauri (Romano 
et al., in preparatio n) considering it as a dwarf galaxy and following the ap- 
proach presented by iRomano. Tosi. Matteucci I ( 2006 . h ereafter, RTM06) fo r 
other dwarfs. We adopt the SF history derived from the ISollima et al] (2005) 



data and allow for both galactic winds and infall of metal poor gas. The model 
predictions are compared with the observational mass, age-metallicity relation, 
metallicity distributions, chemical abundances and abundance ratios. In agree- 
ment with Bekki & Norris results, in no way are we able to obtain model predic- 
tions consistent with the data if we consider uj Centauri as an isolated system. On 
the other hand, by considering it as the residual of a nucleated galaxy stripped 
by the Milky Way 10 Gyr ago, we reproduce rather well all the data, except 
the extremely high He abundance of the blue MS. To achieve this goal, new ad 
hoc assumptions are needed and will be the subject of further efforts, involving 
appropriate assumptions on the fate of the stellar ejecta and new stellar yields 
for both high and intermediate mass stars. Our current, preliminary results are 
shown in FigfTJ 
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3. AGB stars as ISM polluters: 3 He, 12 C and 13 C 

ijj Centauri may require special yields and evolutionary conditions to be ex- 
plained, but the need of improved yields is a much more general problem. Ex- 
amples of important elements for which improved yields are badly needed are 

the He, C, N and O stable isotopes. 

In the late nineties it has been shown (e.g. iGalli et al.1 [l997l : iTosil kOQO. 



and references therein) that to let the predictions of Galactic chemical evolution 
models reproduce the low 3 He abundances measured in Galactic HII regions a 
mechanism is needed, able to drastically reduce the 3 He production normally pre- 
dicted for l ow-mass stars. Such a mechanism was suggest ed by e.g. ICharbonnel 



(1995) and IWasserburg. Boothrovd. fc Sackmann I ( 19951 ) to be the consequence 



of extra-mixing at work in RGB stars, possibly as a consequence of rotation. 
To reconcile the low HII regions abundances with the high 3 He measured in a 
few Planetary Nebulae (PNe), the extra-mixing should affect about 90% of low 
mass stars. Since this extra-mixing implies not only a significant 3 He depletion, 
but also a larger conversion of 12 C into 13 C in such a large fraction of stars, it 
is important to check whether the corresponding yields are consistent with the 
carbon ratios observed in PNe. 



100 



10 - 



/ 



u 



Ejecta 

at AGB tip 



1 2 3 4 5 
Progenitor Mass (M Q ) 



30 



20 

u 

10 



O \ ■ . N7027 

A '\ |Ka:U* | ' A "" , Z = 0.007 

Ml-17 v V'v ' :' - Z = 02 

N8781 A ^f?--"I Z = - 02 

y ~~- z=o.ooi 



1±u 



Surface 
at RGB top 



1 2 3 4 5 
Progenitor Mass (M ) 



Figure 2. 



carbon isotopic ratios measured in PNe 



(dots with error bars in both panels) compared with stellar nucleosynthesis 
predictions. The curves in the left-hand panel s how the r atio p redicted with- 
out deep-mixing just before the PN ej ection bv lMarigol (1200 long-dashed, 
Ivan den Hoek fc Groenewegen I (|1997l ). dotted, and I Forest ini fc CharbonneTI 
(|1997D . dashed. The dotted curves in the right-hand panel show the rati o 
predicted at the end of the RGB phase by iBoothrovd fc Sackman I |l999) , 
with (the curves falling down to low ratios) and without deep-mixing (those 
staying at high ratios). 



In FigElthe carbon isotopic ratios measured by lPalla et al.l ( 20001 ) in several 
PNe are compared with the predictions from various theoretical yields. It is ap- 
parent that the standard nucleosynthesis predictions of the left-hand panels over- 
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produ ce the carbon ratio, while the deep-mixing predictions bv lBoothrovd Sz Sackman" 
(199§) in the right-hand panel nicely fit the data. This nice fit is however mis- 
leading, because Boothroyd & Sackman's computations reach only the RGB 
tip and not the final evolutionary phases. The carbon abundances can be sig- 
nificantly affected by the remaining evolution. Unfortunately, no yields taking 
deep-mixing into account have been calculated beyond the RGB and none of the 
yields computed up to the final phases include the deep-mixing effect. In other 
words, a direct check of the existence of the extra-mixing and of the solution to 
the 3 He problem, with all its cosmological consequences, is currently unfeasible 
! 



4. AGB stars as ISM polluters: nitrogen versus oxygen 

Nitrogen and oxygen are important tracers of the chemical evolution of galaxies. 
Not only are they among the most abundant elements; they are also measurable, 
from the emission lines of HII regions (and PNe in some cases), in galaxies up 
to rather large distances, where they are often the only available metallicity 
indicators. Oxygen is mainly produced by massive stars, whilst nitrogen is 
mostly synthesized by inte rmediate-mass stars. In the framework of the simple 
model of galaxy evolution ( Tinslevlll98dl ) the enrichments of primary elements 



are independent of each other, whilst the growth of a secondary element goes 
as the square of the growth of a primary. The double nature (primary and 
secondary) of the N production beca me apparent when t he N and O measured 
in HII regions of a number of spirals (Diaz &: Tosi 1119860, inc luding the Milky 
Way, and of late- type dwarf galaxies (jMatteucci <fe Tosil Il985l ) were compared 
with the predictions of adequate chemical evolution models: a significant fraction 
(between 30% and 60%) of N must be primary to explain the rather flat trend of 
N/O vs O/H inferred from HII regions in both late-type dwarfs and individual 
spirals. The straight lines in the left-hand pan el of Fig|5] show lea st-squares 
fits to the N/O vs oxygen abundances derived by iDiaz fc Tosi I ([1986ft from the 
HII regions in the Galaxy, M3 1, M33, M101, NGC2403 and IC 342 (for more 
accurate and updated data see iPilvugin. Thuan. fc Vilchez I [2003). The fits for 



these individual spirals are overplotted on a very recent version of the N/O vs 



O/H diagram including results from the Sloan Digital Sky Survey ([Liang et al 
20061 1 



The N/O vs O/H diagram is sometimes interpreted as an evolution plot. 
However, the oxygen abundance in abscissa cannot be taken as a proxy for time. 
In fact, each of the plotted dots corresponds to an HII region, whose abundance 
reflects the final (current) result of the whole history of chemical enrichment in 
the host galaxy region. Such history is known to be quite different for different 
types of galaxies, since it depends on SF history, gas accretion (infall) and loss 
(winds), interactions with companions, etc. As such, the position of each point 
in the plot depends on the different evolutions of different galaxies, and not 
only on the nature (primary or secondary) of the N nucle osynthesis (see also 
Edmunds & Paeel 1 1 19781 : IPilvugin. Thuan. Vilchez II2003T ). 



Much attention has been payed to the metal-poor end of the N/O vs O/H 
diagram, where late-type metal-poor dwarfs are located, to discuss either the 
possibility of primary production of N in massive stars, or the possible connection 
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Figure 3. N/0 vs 0/H in different galaxies as derived fro m HII region 
observ ations (dots). The plot on the left-hand panel is from iLiang et al.l 
(2006) and shows the SDSS dat a as dots. The thick straight lines show the 
trends in individual spirals from lDfaz fc Tosil (|1986f ). The right-hand panels 
are blow-ups of the metal-poor region of the diagram and the lines show some 
of the effects of galaxy evolution (see text for details) . 



between late- type dwarfs and Damped Lyman-a systems. The plots on the 
right-hand side of Figj3] show some of the effects on the N and O abundances 
occurring during galaxy evolution. The lines shown in the two diagrams on the 
right correspond to the predictions of different chemical evolution models for 
the time behaviour of oxygen and nitrogen in a late-type dwarf galaxy. The 
top-right panel compares with the abundances measured fro m the HII r egion s 
observed in late-type dwarfs (dots) the model predictions by iPilvugin I ( 19931 ) 
for a generic Blue Compact galaxy experiencing about ten SF bursts. It shows 
how the HII region self-enrichment and the different timescales for the stellar 
ejection and subsequent diffusion in the ISM of oxygen and nitrogen lead to a 
saw-tooth shape for the N/O ratio as a function of oxygen. 

The bottom-right panel show instead the predictions by RTM06 for the 
starburst dwarf NGC 1569 (always the line with higher N/O for each line-type; 
blue in the colour version of the figure) and the Blue Compact dwarf NGC 1705 
(always the line with lower N/O for each line- type; red in the colour version of the 
figure). The big symbols with error bar indicate the corresponding values derived 
from HII regions in NGC 1569 (blue square) and NGC 1705 (red diamond). 
The other points correspond to HII regions values in other late-type dwarfs. 
Different line-types correspond to different assumptions on the stellar yields or 
on the galactic wind efficiency. See RTM06 for references and details. All the 
lines show model predictions for the evolution of oxygen and nitrogen with time 
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and it is apparent how much the time behaviour can be different from a simple 
fit to the current N/O vs O/H distribution inferred from HII regions. 

In spite of the different clas sification, NGC 1569 and NGC 1705 have 
quite similar properties (see e.g., lAnnibali et al.l [20031 ; lAngeretti et al.l 2005. 



and references therein): similar gas and star masses, metallicity, IMF, SF his- 
tory with very strong recent activity, and observational evidence of similarly 
strong galactic winds. Yet, they have quite different N/O ratios, which re- 
quire different evolutionary assumptions. The long-dashed lines in the bottom- 
rig ht panel of FigfS] show the model pre dictions based on the standard yields 
by Ivan den Hoek Sz Groenewegen" (1997) for intermediate mass stars. Clearly 
a standard N production may be consistent with the N/O ratio observed in 
NGC 1569, but un-reconcilable with the low N/O of NGC 1705. If we con- 
sider yields with smaller N production, even in the extreme case of those by 
iMevnet fe Maeder ( 20021 ) where the hot-bottom burning phase has not been 
computed, the models overpredict the N/O ratio observed in NGC 1705 (short- 
dashed lines). A better fit to the observed abundances is obtained (RTM06) if 
not only the nitrogen production in intermediate mass stars is relatively low, 
as predicted for instance by the minim al hot-bottom burning case proposed by 
van den Hoek fe Groenewegen I ( 19971 ). but also a higher efficiency of nitrogen 



loss in the galactic wind is allowed in NGC 1705 than in NGC 1569 (solid lines) 
Detailed chemical evolution models of individual dwarfs, such as those 
shown here for NGC 1569 and NGC 1705, have become possible only recently, 
mostly thanks to HST, which has allowed to derive their SF histories back 
to quite early epochs. Two kinds of models for individual galaxies can be 



comp uted: standard chemical evolut ion models fe.g. ICarigi. Colin. Peimbert . 
ll999l : lLanfranchi fc Matteucci ll2003|:lRomano. Tosi. Matteucci 1 12006) and chemo- 



dynamical models fe.g iRecchi et ai1l200l l2006 ) . The former have the drawback 



of a simplistic treatment of star and SN feedbacks and gas motions, the latter 
have the problem that the timescales appropriate for hydrodynamics make it 
prohibitive, in terms of CPU time, to follow the system evolution over more 
than 1 Gyr. The challenge in the next few years is to improve both types of ap- 
proaches and get a more realistic insight of how stars and gas evolve, chemically 
and dynamically, in their host galaxies. 



5. The need for improved yields 

In the previous sections, some examples of the uncertainties affecting the pre- 
dictions of chemical evolution models have been presented. Part of these uncer- 
tainties are due to the lack of complete and homogenous sets of stellar yields 
for various initial metallicities. To compute adequate chemical evolution models 
of whatever galaxy, we need homogeneous chemical yields for all the major iso- 
topes, for the whole range of stellar masses, for many initial compositions and 
taking into account all the most relevant processes occurring in the stellar inte- 
riors until the final evolutionary phases. Such optimal grid of yields is far from 
existing. The results from stellar nucleosynthesis are steadily improving with 
time, with most of the processes occurring during stellar evolution being treated 
with increasing precision. However, the yields available to the community are 
still very heterogeneous and incomplete and not one single set of nucleosynthesis 



8 



Tosi 



predictions exists taking properly into account all the processes in all the phases 
of stars of all masses (say from 0.8 to 100 M ) and at least 2-3 different metal- 
licities (from metal poor to solar and, possibly, super solar). This circumstance 
not only prevents the computation of detailed self-consistent chemical evolution 
models for a number of key elements, but can even lead to misleading results. 
The potential risk can be visualised by comparing with each other the yields 
provided by different authors and the corresponding normalizations through an 
IMF. 




Figure 4. Stellar yields: ejected mass of newly synthesized element as a 
function of the stellar initi al mass. Left-ha nd p anel: the case of the homoge- 
neous sets of solar yields bv lMarigol (|200ll ) and lPortinari. Chiosi. fe Bressan I 
(1998) covering the whole mass range. Right-hand pan el: the stan- 
dard case of incomplete c overage with inhomogeneous sets (|Marigo 1 120011 : 
IWooslev fc WeaverTll995f ) which do not consider initial masses with 5 < 
M/M Q < 11 (the lines plotted in this mass range are linear interpolations) 
and beyond 40 M Q . 



The left-hand panel of Fig j] show s a subset of the solar yields presented 
bv lPortinari. Chiosi. &: Bressan ( 19981 ) for massive and quasi-massive stars and 
by iMarigo (2001) for low and intermediate mass stars. This is the best case 
in literature of self-consistent yields for all masses, based on the same stellar 
evolution models and input physics. What is normally available in the literature 
is a collection of yields for partial mass ranges, each computed under different 
assumptions and often for different metallicities. A typical case is shown in the 
right-hand panel of FigU where the solar yields bylMjmgo" ( 2001) for low and 
i nterm ediate mass stars are now combined with those by IWooslev &: Weaver 



(|1995T ) for massive stars. Two problems are immediately apparent. First, stars 
in the mass range 5 < M/Mq < 11 have not been computed by either Marigo 
or Woosley & Weaver, which implies that to use this combination of yields 
one must interpolate over this mass interval. Second, massive stars do not 
go beyond 40 M Q and to higher masses one must therefore extrapolate. The 
latter issue might not have overwhelming consequences in the modelling of the 
recent chemical evolution of the Milky Way, since any reasonable IMF predicts 
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Figu re 5. Fractio nal ejected mass of newly synthesized element weighted 
with lTinslev I (1980) I MF. L eft-ha nd p anel: the case of the homogene o us set s 
of solar yields by iMarigo I (|2001[) and iPortinari. Chiosi. fc Bressan I (| 19981) . 
Right-hand pa nel: the case of the inhomog eneous, non adjacent sets by 
IMarigo I (pOOll ) and lWooslev fc Weaver I (|1995l ). 



very few stars more massive than 40 Mq, but can be extremely relevant for 
very early epochs, when the most massive stars were the only polluters. The 
former problem has very serious implications because stars in the 5-11 M Q 
range are the most effective contributors to the ISM chemical enrichment. In 
FigGD the yields of FigH have been weighted with Tinsley's (1980) IMF. The 
linear interpolation performed to cover the 5-11 M interval absent in the 
Marigo/Woosley& Weaver combination results in a bump (right-hand panel) in 
the contribution of these stars to the enrichment of He, N and O which is totally 
absent in the left-hand panel, where the homogeneous sets of yields are shown. 
This enhanced contribution is most likely spurious and can lead to a significant 
overprediction of the elements mostly produced by stars of these masses. 

To overcome these problems, we strongly encourage the community of stellar 
nucleosynthesis experts to provide homogeneous yields for all stellar masses, 
computed up to the final evolutionary phases and for several initial metallicities. 
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